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ABSTRACT 

This is a part of a multiwavelength study aimed at use of complementary photometric, 
polarimetric and spectroscopic data to achieve an understanding of the activity process 
in late-type stars. Here we present the study of FR Cnc, a young, active and spotted 
star. 

We performed analysis of AS AS — 3 (The All Sky Automated Survey) data for the 
years 2002-2008 and amended the value of the rotational period to be 0.826518 d. The 
amplitude of photometric variations decreased abruptly in the year 2005, while the 
mean brightness remained the same, which was interpreted as a quick redistribution of 
spots. BVRc and Ic broad band photometric calibration was performed for 166 stars 
in FR Cnc vicinity. 

The photometry at Terskol Observatory shows two brightening episodes, one of 
which occurred at the same phase as the flare of 2006 November 23. Polarimetric 
BVR observations indicate the probable presence of a supplementary source of po- 
larization. We monitored FR Cnc spectroscopically during the years 2004-2008. We 
concluded that the RV changes cannot be explained by the binary nature of FR Cnc. 
We determined the spectral type of FR Cnc as K7V. Calculated galactic space- velocity 
components {U, V, W) indicate that FR Cnc belongs to the young disc population and 
might also belong to the IC 2391 moving group. Based on Li i A6707.8 measurement, 
we estimated the age of FR Cnc to be between 10-120 Myr. Doppler Tomography 
was applied to create a starspot image of FR Cnc. We optimized the goodness of fit 
to the deconvolved profiles for axial inclination, equivalent width and v sin i, finding 
V sin i = 46.2 km s~^ and i — 55°. We also generated a synthetic V-band lightcurve 
based on Doppler imaging that makes simultaneous use of spectroscopic and photo- 
metric data. This synthetic lightcurve displays the same morphology and amplitude 
as the observed one. 

The starspot distribution of FR Cnc is also of interest since it is one of the latest 
spectral types to have been imaged. No polar spot was detected on FR Cnc. 

Key words: stars: activity - stars: flare - stars: rotation - stars: individual: FR Cnc. 
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1 INTRODUCTION 

It is well-known that late-type stars show magnetic activity 
similar to the activity of our Sun, but the physics of 'stellar' 
activity is not yet well understood. In addition, the activity 
level manifested by late-type stars is much higher than that 
observed for the Sun. 

FR Cnc (= BD-H16 1753 = MCC 527 = lES 0829+15.9 
= IRXS J083230. 9-^154940 = HIP 41889) was first men- 
tioned as a probable active star when it was identi- 
fied as the optical counterpart of a soft X-ray source 
lES 0829-1-15.9 in the Emstem Slew Survey. It has 
V = 10.43 mag, spec tral type K8V, the X-ray flux of 
«10"" erg s"^ cm"^ (|Elvis et al.lll992l : ISchachter et all 
1 19961 ). Lately, this object was rediscovered as an X-ray 
source IRXS J083230.9-f 154940 in the ROSAT All-Sky 
Survey (RASS) with lowe r X-ray flux at the level of 

2 X 10"^^ erg s"^ cm"^ (|Voges et all 1 19991 ). The X-ray 
luminosity of (2 — 12) x 10^^ erg s~^ and the 
tio of X-ray to bolometric luminosity of ^ 10 
(|Pandev et al.l |2005|) indicates th at this object has an ac- 
tive corona ( Schachter et al.lll996l ). 

In the Hipparcos catalogue this star was mentioned as 
an unsolved variable star with the i dentifier HIP 41889 an d 
0.17 mag amplitude of variability (|Perrvman et al.lll997^ . 
It was classified as BY Dra type star (i.e. its variability 
caused by rotational modulation of sta r spots) and given 
the name FR Cnc by iKazarovets et all (1 19991 ). For anal 
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ysis of Hipparcos observations see iPandev et all (|2005h . 
FR Cnc (RA(2000) = 08''32™3O!5287 and Dec. (2000) = 
+ 15°49'26"193) has 30.24 ± 2.03 marcsec parallax 
jPerrvman et al.lll997l ) that implies a distance of 33 ± 2 pc 
and an absolute magnitude of 7.8. The kinematics of FR Cnc 
suggests that it is a very young (35-55 Myr) main-sequence 
star and a possi ble member of the IC 2391 supercluster, as 
it was shown by IPandev et all (j2005h . lUpgren et all (j2002l ) 
concluded that FR Cnc is not a binary system, based on two 
measurements of the RV. 

The presence of Ca ll H & K and Ha emission lines in 
the spectra indicates high chromos pheric activity in FR Cnc 
(jPandev et al.ll2002l : IPandevllioO^ ) . In 'quiescent' state this 
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Figure 1. Lightcur ve of FR Cnc on 2006 November 23 from 
iGolovin et all | |2007|) : " The flare of FR Cnc: shifted differential 
lightcurves in B—, V—, R— and /—bands as well as the difi'erence 
check star — comparison star ('Ch' on the plot)" 



object manifests optical variability with the dominant pe- 
riod 0.8267 ± 0.0 004 d due to the p resence of starspots 
and axial rotation (|Pandev et al.l [20051 ). In addition, pho- 
tometry obtained in 2005 February - April with Kilodegree 
Extremely Little Telescope indicates FR Cnc optical vari- 
ability with 0.827 d period (jPepper et al.|[2008l ) when mon- 
itoring the Praesepe open cluster for transiting exoplanets. 

The first ever-detected optical fiare of this object 
was observed during CCD photometry of FR Cnc on 
2006 November 23 at Crimean Astrophysical Observatory 
(U kraine) with 38-cm Cassegrain telescope and described 
bv lColovin et all (|2007l ). The fiare was observed in BVRI- 
bands (see Fig. [T|. The amplitude reached even 1 mag in the 
_B-band and was decreasing towards the /-band. The fiare 
energy output in the _B-band was about 1.73 x lO'^^erg A^^ 
and fiare to quiescent fiux ratio was 38.63 per cent. 



2 PHOTOMETRIC OBSERVATIONS 

Most of the information on the photospheric activity (e.g. 
starspots) of BY Dra-type stars comes from photometric 
observations. The mean brightness level is strongly depen- 
dent on the percentage of spotted area of the surface, while 
changes in spot distribution over the surface could result in 
changes of the amplitude of variability. FR Cnc has a short 
(for such class of objects) rotational peri o d of .8267 d. As 
was shown in iDorren. Guinan fc Dewarj (| 19941 ) , this short 
rotational period leads us to expect large fiare activity of 
the star. 

The detection of a fiare on FR Cnc on 2006 November 
23 motivated us to continue photometric monitoring of this 
object as well as to study its archival ASAS-3 observat ions 
(The All Sky Automated Survey; see lPoimans 3 (|2002l ) for 
description of equipment and data pipeline). 



2.1 AS AS Photometrical Observations 

FR Cnc was observed in l/-band with ASAS-3 survey dur- 
ing 2002 December - 2008 May (6 observational seasons, see 
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Figure 2. ASAS long term light curve 

Fig. ID). AH the data were split into separate data sets accord- 
ing to the 'seasonal ga ps' in observations an d folded with the 

0. 8267 d period ffrom IPandev et al.llioosi ') and represented 
in Fig. [3] (plotted twice for clarity). Table [T] represents the 
log of observations. 

The initial epoch was common to calculate phases 
for all 6 phase diagrams and was chosen arbitrary as 
HJD (UTC) = 2452635.72669 (first point in dataset). No 
evidence of flares in the ASAS-3 data was found. The verti- 
cal dashed line on Fig. [3] indicates the phase when the flare 
on 2006 November 23 occurred. 

ASAS-3 data are not covering flare on 2006 November 
23. Nearest ASAS-3 observations were done on JD 2453881 
and JD 24544091, what is 181 days before and 29 days af- 
ter the flare. ASAS-3 observations in 2006-2007 at a phase 
~0.88 (143 days after the flare) show that FR Cnc was 
brighter than during the rest of the time. This probably 
could be related to brightening episode, which was detected 
at the same phase during Terskol observations (see sect. 12. 3( ). 

Profiles of variability as well as the amplitude of vari- 
ability are different from season to season, while the mean 
brightness remains constant (Vmcan = 10.439 mag within the 
error limits of ct = 0.017 mag). To illustrate this, we plot- 
ted the amplitude and the mean brightness as a function of 
mean epoch of observations (Fig. |3| . An abrupt decrease of 
amplitude in the year 2005 is clearly seen. One of the possi- 
ble interpretations is continuous spottedness of the star and 
redistribution of spots/spot groups from season to season: 

1. e. spots didn't disappear, but distributed more uniformly 
over FR Cnc surface. To support this idea, notice the con- 
stancy of the mean brightness level. If spots disappear then 
increase of brightness and decrease of amplitude could be 
expected. 



2.1.1 Penodogram Analysis 

The data set was searched for periodic variation of bright- 
ness (in order to estimate with better accuracy the known ro- 
tational period) using the period04 package, developed by 
Patrick Lenz (Institute of Astrono my, University of Vienna, 
Austria; see iLenz fc BregeJl2005l ). Discrete Fourier trans- 
form (DFT) algorithm was applied for statistical analysis. 
Julian dates are heliocentrically corrected. The average zero 
point of 10.439 mag was subtracted to prevent the appear- 



ance of additional features on the periodogram centered at 
frequency 0.0. 

Making a periodicity analysis of ASAS-3 photometry, 
we found a dominant frequency of / — 1.209895 c d^^, 
while, as defined by equation [TJ cr(/) = 0.000022 (there- 
fore, period P = 0.826518 ± 0.000015 d; see Fig. [S]). The 
detected periodicity could be interpreted as the rotational 
pe riod and it is i n good agreement with the period founded 
bv iPandev et al.] (|2005l ). but with improved accuracy (due 
to the longer time span of ASAS-3 observations). No other 
periodicity of FR Cnc brightness modulations was found on 
the basis of ASAS-3 observations. 

The DFT routine was applied separately to each season 
of observations as well. The obtained periodograms do not 
reveal any other significant periodicity, but only the same 
peak as for the periodogram for the whole time-string. The 
dominant frequency remains constant within the error limits 
of o"(/) = 0.0003 for all of 6 periodograms. 



2.1.2 Estimation of Accuracy and Reliability of the 
Detected Period 



Emp irical results from observati onal analvsis (|Breger et al.l 
Il993h and numerical simulations (|Kuschnig et al.lll997 ) have 
shown that the ratio in amplitude in the periodogram be- 
tween signal and noise should not be lower than 4.0 to give 
good confidence in the detected peak. We calculated the 
S/N-ratio from the periodogram for the determined domi- 
nant frequency, so S/N = 6.88. 

The Spectral Window Function is important to confirm 
that the obtained frequencies are real or an artefact of the 
window function. This was done by assigning 1 to bright- 
ness values with the same observation times and checking 
the resulting diagram. There was no evidence of significant 
power at the location of the peak (Fig. [5]). The dominant 
frequency on the spectral window is _F = 1.0027 with ampli- 
tude A — 0.9391, which is due to daily gaps in observations. 

Parameter uncertainties were calculated from an error 
matrix, which is a by-product of non-linear least-squares fit- 
ting procedure. Other types of uncertainties are those which 
could be calculated from analytically derived formulae as- 
suming an ideal case. Based on some assumptions one can 
derive a formula for the uncertai nties in frequency and signal 
ampli tude at this frequency. Se e lMontgomerv fc Odonoghue 
( jl999[) and iBreger et al.1 (|l999l 'l. iKallinger. Reegen fc Weiss 
(2003) for thfi derivation based on a monoperiodic fit. The 
determined equation can be applied for each frequency sep- 
arately; 




(1) 



(2) 



where A*' is the number of time points, T is the time 
length of the data set, a{m) denotes the residuals from 
the fit and a refers to the signal amplitude at the particu- 
lar frequency. Both 'analytical' and 'least-squares error ma- 
trix' calculations give similar results: cr(/) — 0.000022 and 
a{a) = 0.0037. So, on the periodogram (Fig. [S} we plotted 
the 3(T(a)-level (dashed line) to show that the detected peak 
exceeds it significantly. 
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Figure 3. ASAS-3 phase diagrams for six observing seasons. Please note that scales are the same for all plots. Vertical dashed line 
indicates the phase where the flare on 2006 November 23 was detected (see text for the explanation). 



Another method for Icr-level of a mplitude calculation 
was p roposed for use by R.A. Fisher ( see lFisheiiri935l : iFisheJ 
Il936h and often called as Fisher Randomization Test. The 
idea is to take the original hght curve and preserving the 
time column, shuffle the corresponding intensities around. 
That destroys any coherent signal in the hght curve while 
keeping the time sampling intact. We are left with a shuffled 
light curve of pure white noise. The next step is to com- 
pute a DFT of this light curve which will look really noisy. 
The standard deviation of the average amplitude of such 
a DFT is close to the la limit. We iterated this 25 times, 
therefore we use the average of 25 standard deviation val- 
ues to determine la-level. This pipeline gave us the value of 
o"(a) = 0.0028, which is slightly less than a-{a) from 'analyt- 
ical' and 'least-squares error matrix' calculations, hence we 
plotted the bigger value on our periodogram to be confident. 
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Figure 4. Mean amplitude of variability, EW of Ha and mean 
brightness of PR Cnc during the years 2002-2008. (See sect. [6711 
regarding Ho.) 
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Figure 5. Periodogram of ASAS-3 data in frequency range of 1.0-1.5 and Spectral Window 



Table 1. Log of ASA5-5-obscrvations of FR Cnc 



Year 


Tstart 


Tend 






2450000+ 


2450000+ 




2002-2003 


2635 


2791 


59 


2003-2004 


2924 


3146 


62 


2004-2005 


3357 


3512 


63 


2005-2006 


3674 


3881 


61 


2006-2007 


4091 


4247 


34 


2007-2008 


4409 


4611 


56 



Table 2. Log of Tcrskol photometric observations of FR Cnc in 
the years 2007-2008 



2.2 Photometric Sequence 

We carried out the BVRcIc photometric cahbration for 166 
stars in the vicinity of FR Cnc with l/-magnitudes in the 
range 9.85-18.06 mag that could serve as comparison stars. 
Calibration was done at the Sonoita Research Observatory 
(Arizona, USA) using a robotic 0.35-m telescope, equipped 
with an SBIG STL-IOOIXE CCD camera. A table with this 
photometrical sequence is available electronically only via 
the AAVSO ftp-serveiQ. For user convenience , we have used 
the ALADIN Sky Atlas (|Bonnarel et al.l [2OO0I ') to align our 
calibration data on a DSS2/STScI POSSII imagfEI. 



2.3 Terskol Observations 

Optical B-band photometry was carried out from 2007 
March to 2008 February at Terskol Branch of the Astronomy 
Institute (Russia) with 29-cm telescope and Apogee-47 Alta 
CCD camera. All observations were made in the B-band as 
the flare amplitude is expected to increase with decreasing 
wavelength. The duration of each observing run varies from 
2 to 7 h. See Table [2]for log of observations. The calibration 
process of the obtained frames, c omparison and check stars 
remains the same as described bv lColovin et al] (|2007t ). 

The mean amplitude of FR Cnc brightness variations 
in the B-band during the observations in the year 2007 was 
0.12 mag and 0.13 mag in the year 2008. There are no pecu- 
liar features in the light curve from the year 2008, while in 
the photometry obtained in the year 2007 two brightening 

^ ftp:/ /ftp. aavso.org/public/calib/frcnc.dat 

^ http:/ /www. mao.kicv.ua/ardb/ref/agolovin. html 



Year 


Beginning of the Rliu 
(H,ID) 


End of the Run 
(H,ID) 


2007 March 


2454171.2453 


2454171.4361 




2454174.3501 


2454174.4019 




2454180.1863 


2454180.4691 




2454182.2395 


2454182.4972 




2454188.4424 


2454188.4721 


2008 February 


2454498.4150 
2454500.2675 
2454501.1575 
2454502.1681 
2454503.2157 
2454504.1755 


2454498.5233 
2454500.5345 
2454501.5403 
2454502.5870 
2454503.4481 
2454504.3025 



episodes were detected: namely, on 2454180.3 and 2454182.3 
(see Fig. |6}. Amplitude in B-band was 0.06 and 0.12 mag re- 
spectively. It has to be noted that the second episode (HJD 
= 2454182.3) occurred at the same phase as the flare of 2006 
November 23 (phase = 0.88). Probably both the events (flare 
on 2006 November 23 and brightening episode at phase 0.88) 
are originated from the same long-living active regions on 
the surface of FR Cnc. 



3 POLARIMETRIC OBSERVATIONS 

The BVR broad-band polarimetric observations of FR Cnc 
were obtained on 2007 October 19 and 20 using TK 1024 
pixel^ CCD camera mounted on the Cassegrain focus of the 
104-cm Sampurnanand Telescope of ARIES, Nainital (In- 
dia). The optical imaging polarimetry was carried out in 
B, V and R [Xs^a = 0.440 /im, Av^j.,. = 0.550 /im and \r^„ 
= 0.660 ^m,) photometric bands. 

Standard stars for null polarization and for the zero- 
point of the polariza tion position angle were taken from 
ISchmidt et all (|l992t ). The results for standards are given 
in Table [3] From the results, it can be concluded that ob- 
tained values of polarization and pos ition angles are in good 
agreement with lSchmidt et al.l (|l992l ) within the error limit. 

Both the program and the standard stars were observed 
during the same night. HD 25433 and HD 19820 (= CC Cas) 



6 A. Golovin et al. 




I I I I I I I 

0,0 0,5 1,0 1,5 2,0 

Phase 

Figure 6. _B-band Photometry at Terskol Observatory in 2007 Marcli. Plotted twice for clarity. Dotted lines denote phases when 
polarimctric observations were done. Note an enlarged plots of brightening episodes on JD = 2454180 and JD = 2454182. 



Table 3. Observed polarized standard stars. 



Star 


name 


Filter 


Piep {percent) 




Pitep (percent) 










Published data 


This 


paper 


HD 


25433 


B 


5.23±0.09 


134.3±0.05 


5.17±0.21 


135.6 ±1.0 






V 


5.12±0.06 


134.2±0.03 


5.13±0.09 


133.5 ±0.8 






R 


4.73±0.05 


133.6±0.03 


4.76±0.13 


132.9 ±0.5 


HD 


19820 


B 


4.70±0.04 


115.70±0.22 


4.66±0.07 


115.49 ±0.19 






V 


4.79±0.03 


114.93±0.17 


4.76±0.10 


114.15 ±0.20 






R 


4.53±0.03 


114.46±0.17 


4.56±0.17 


114.18 ±0.21 



were used as a standard polarized stars, while HD 21447 and 
G 191-B2B (= HIP 23692) served as standard unpolarized 
stars. 

The results are listed in Tables [3l |4] and [5l where P is 
the fraction of the total hght in the linearly polarized condi- 
tion and 6 is the position angle of polarization plane to the 
equatorial plane. It is denoted by the normalized Stokes' pa- 
rameter q (= Ql I), when the half wave plate's fast axis is 
aligned to the reference axis (a = 0°). Similarly, the nor- 
malized Stokes' parameter u {= U/I), when the half wave 
plate is at 22.5°. For further details on used equipment and 
the method of observations, refer to lMedhi et al.l ( )2007l l. 

By dotted lines on the lightcurve we indicate the phases 
(V3i = 0.16, V32 = 0.13 for 2007 October 19 and 20 respec- 
tively), when our polarimetric observations were conducted 
(see Fig. El . FR One was in maximal brightness during that 
time. 

Polarization in FR Cnc could be magnetic in ori- 
gin. The degree of polarization depends nonlinearly on 
the size of magnetic regions (see iHuovelin fc Saaij Il99ll : 
ISaar fc Huovelinill993l ). These authors have also calculated 
a grid of expected degrees of polarization in UBVRI band 
for stars with temperature from 4000 to 7000 K and log g 
from 2.0 to 4.5. We have used their results to compare our 
observed values of polarization for FR Cnc. Fig. [8] repre- 
sents the degree of polarization for FR Cnc in BVR bands. 



The maximum possible degree of polarization for the total 
spot area of 24 per cen t is derived from the calculations of 
ISaar fc HuovelinI (jl993l ) for the star corresponding to the 
spectral type of K5-7V and characteristic magnetic field of 
2.7 kC These values are over-plotted and represented by a 
solid line in Fig. [HI The observed polarization in the B band 
is in good agreement with the theoretical values expected 
for Zeeman polarization model. However, the observed po- 
larization in V and R bands slightly exceeds the theoretical 
values. 

Model values for K2V-K7V and K2IV-K2V spectral 
types appear to be even lower and certainly do not match 
the polarization in any of observed bands. This is also ob- 
served in so me other young spotted stars (MS Ser, LQ Hya, 
VY Ari; see lAlekseevll200a ) and probably due to the pres- 
ence of a supplementary source of linear polarization. 

The predicted values of polarization due to Thompson 
and Rayleigh scattering from inhomogeneous regions are not 
enough to explain the observed polarization excess (Thomp- 
son and Rayleigh scattering for the assumed spectral type 
suppos ed to be of order of 10 ~^ and 10~* per cent respec- 
tively; ISaar fc HuovelinI 1 19931 '). The mechanism which can 
produce additional linear polarization is probably scatter - 
ing in circumstellar material (e.g. see iPandev et al]|2009l ): 
on the other hand, the mentioned models are unacceptable if 
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Figure 7. B-band Photometry at Tcrskol Observatory in 2008 
February. 



Table 4. Observed unpolarized standard stars 



Star name 


Filter 


<3 


u 


HD 21447 


B 


0.019 


0.011 




V 


0.037 


-0.031 




R 


-0.035 


-0.039 


G191B2B 


B 


0.072 


-0.059 




V 


-0.022 


-0.041 




R 


-0.036 


0.027 



FR C nc is a close binary star |Alekseevl20o3 : lElias fc DorrenI 
Il99d and[s aar fc Huovelinlll993h . 



4 SPECTROSCOPIC OBSERVATIONS 

A total of 58 high and low resolution spectra of FR Cnc have 
been obtained and analysed in this work. The spectroscopic 
data were obtained during five observing runs. Details of 
each observing run are given in Table [G] date, telescope, 
spectrograph, CCD chip, spectral range covered, number of 



c 
u 

a. 

0. 



1.2 



0.9 



0.3 



0.0 



A 19 October 2007 
O 20 October 2007 
■ Model values for K5-7V stars 



0.4 



0.5 

X (iivn) 



0.6 



0.7 



Figure 8. The degree of polarization of FR Cnc as a function 
of waveleng th. Model values of po larization for K5-7V spectral 
types stars llSaar fc Huovelinl . ll993fl are plotted for comparison. 



orders included in each echelle spectrum, range of reciprocal 
dispersion, spectral resolution (determined as the full width 
at half maximum, FWHM, of the arc comparison lines) and 
mean S/N in the Ha line region. 

The spectra were extracted using the standard reduc- 
tion procedures in the iRAljfl echelle package (bias subtrac- 
tion, flat-field division and optimal extraction of the spec- 
tra). We obtained the wavelength calibration by taking spec- 
tra of a Th-Ar lamp. Finally, we normalized the spectra by 
a polynomial fit to the observed continuum. 



5 STELLAR PARAMETERS 

Stellar parameters of FR Cnc are given in Table [7| and Ta- 
ble [3 The photometric data {B - V, V), -Pphot, projected 
rotational velocity (vsini), and galactic space- velocity com- 
ponents ([/, V, W) have been determined in this paper. The 
astrometric data (parallax, tt; proper motions, UgCOs S an d 
fig) are from Hipparcos and Tycho-2 catalogues (|ESA||1997| ). 



5.1 Spectral classification 

FR Cnc is classified as a K8V-star by lSchacht er ct al.l (ll996l ). 
while a multicolour photometric study allowed Pandev et al.l 
(|2005h to classify it as a K5V-star. They also have obtained 
the Spectral Energy Distribution (SED) of the star matching 
a Teff of 4250±250 K and a log g of 4.50±0.5, that agrees 
with a K5V classification. 

We have compared our high resolution echelle spec- 
tra, in several spectral orders free of lines sensitive to chro- 
mospheric activity, with spectra of inactive reference stars 
of different spectral types and luminosity classes, observed 
during the same observing run. This analysis makes use 
of a modified version of the program STA RMOD (jstar- 
mod) developed at Penn State University jBardenl Il985l : 
iLopez-Santiago et al.|[201ol ) . This program constructs a syn- 
thesized stellar spectrum from artificially rotationally broad- 
ened, radial- velocity shifted, and weighted spectra of appro- 
priate reference stars. For FR Cnc, we have obtained the 
best fit with a K7V reference star, which is in closer agree- 
ment with the K8V classification rather with the K5V from 
iPandey et all l|2005h . 

5.2 Radial Velocity 

We have determined the heliocentric radial velocities (RV 
her eafter) by making u se of cross-correlation technique (see 
e.g. lCalvez et al.]|2007l ). The spectra of the target were cross- 
correlated order by order, using the routine FXCOR in iraf, 
against sp ectra of RV standard s with similar spectral type 
taken from lBeavers et al.l(|l979l ). We derived the RV for each 
order from the position of the peak of the cross-correlation 
function (CCF) and calculated the uncertainties based on 
the fitted peak height and the antisymmetric noise as de- 
scribed by lTonrv fc Davis! l|l979t ). 

^ IRAF is distributed by the National Optical Observatory, which 
is operated by the Association of Universities for Research in 
Astronomy, Inc., under contract with the National Science Foun- 
dation. 
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Table 5. Observed BVR polarization values for FR Cnc. 



Date of Observation 


Filter 


Time (UT) 


PdieiJ (percent) 


e±Efl(°) 


October 19, 2007 


B 


22:22:12.0 


0.22 ±0.05 


57 ±7 




V 


22:01:45.1 


0.20 ±0.02 


55 ±2 




R 


22:10:52.6 


0.16 ±0.05 


61 ±7 


October 20, 2007 


B 


21:25:36.6 


0.26 ±0.04 


55 ±5 




V 


21:30:00.5 


0.24 ±0.07 


54 ±6 




R 


21:35:48.6 


0.17 ±0.05 


58 ±7 



Table 6. Description of spectroscopic observations 



Number 


Date 

(dd/mm/yyyy) 


Telescope 


Instrument 


Detector 




Sped 

(A) 


. range 


Orders 


Dispersion 

(A) 


FWHM^ 
(A/pixel) 


S/N 
Ha 


1 


29/03-07/04/2004 


2.2-m" 


FOCES'' 


2048x2048 24^ SITE#ld 




3450- 


10700 


112 


0.04-0.13 


0.08-0.035 


40 


2 


11-13/04/2004 


2.0-m'' 


HFOSC 


2000x4000 SiTe ST-002, grism 


Grl4 


3270- 


6160 


1 


3.57 


7.23-1.47 


200 










2000x4000 SiTe ST-002, grims 


Gr8 


5800- 


8350 


1 


3.23 


3.66-3.99 


200 


3 


16-21/12/2006 


2.2-m" 


FOCES'' 


2048x2048 24^ SITE#ld 




3600- 


10700 


106 


0.08-0.1 


0.08-0.04 


40 


4 


24-26/02/2007 


2.2-m" 


FOCES'' 


2048x2048 24^ SITE#ld 




3600- 


10700 


106 


0.04-0.13 


0.07-0.41 


40 




07-08/05/2007 






















5 


21/03/2008 


NOT'^ 


FIES^' 


2000x2000 EEV42-40 




3620- 


7360 


80 


0.02-0.04 


0.05-0.11 


80 



^ Full Width at Half Maximum of the arc comparison lines; 

" 2.2-m telescope at the German-Spanish Astronomical Observatory (CAHA, Almeria, Spain); 

^ 2.0-m Himalayan Chandra Telescope at the Indian Astronomical Observatory (Mt. Saraswrati, Hanle, India); 

Nordic Optical Telescope (NOT) at the Observatorio del Roque de los Muchachos (La Palma, Spain); 

The Fibre Optics Cassegrain Echelle Spectrograph (FOCES); 
" Himalaya Faint Object Spectrograph and Camera (HFOSC); 
f The high-resolution Fibre- fed Echelle Spectrograph (FIES). 



As Fig. [9] (top) shows, the irregular profiles of the CCF 
(double peaks and asymmetries) can produce significant er- 
rors in RV measures. These features show regular variations: 
a double peak moving on time-scale of the rotational period 
can be seen. Photospheric activity features on the stellar sur- 
face that disturb the profile of the photospheric lines could 
induce variations in the peak of the CCF, but a stellar com- 
panion could also produce the double peak effect. 

We checked if these RV variations could be due to a 
binary nature of FR Cnc. We found no evidence for the 
existence of a companion by measuring the RV of the peaks 
and by trying to fit the data to a coherent orbit. 

We also carried out a line bisector analysis to enable us 
to ascertain whether the RV variations may be attributed 
to starspots. 

The CCF was computed for regions which include the 
photospheric lines commonly used in the Doppler Imag- 
ing te chnique, while exc luding chromospheric and telluric 
lines jOueloz et al-llioOll ). We computed the bisector and, 
to quantify the changes in the CCF bisector shape, also 
the bisector inverse slope (BIS). The BIS was defined as 
the difference of the average values of the top and bot- 
tom zones (we avoided wings and core of the CCF profile, 
due to errors of bisectors measurements, which are large in 
these zones) . We studied the bisector variations only for FO- 
CES04 run, as it was the more suitable data for the study. 
Th ree cases are repor t ed in bibliography (see for exam - 
ple iQueloz et al.l (l200ll ). lMartfnez Fiorenzano et all (|2005h . 
iBonfils et all (|2007l ). etc.): a) anticorrelation, which indi- 
cates that the RV variations are due to stellar activity (by 
active regions at the stellar surface like spots or plages), b) 



lack of correlation, which indicates the Doppler reflex mo- 
tion around the center of mass due to other bodies orbiting 
the star, c) correlation, which, as pointed out by Martinez- 
Fiorenzano et al. (2005), indicates that the RV variations 
are due to fight contamination from an unseen stellar com- 
panion. 

As shown in Fig. IIOI there is an anticorrelation be- 
tween BIS and RV, with a Pearson correlation coefficient 
(r) of -0.6851. This result suggests that the RV variations 
of FR Cnc are due to stellar activity variations (e.g. spots 
on photosphere) and not due to a binary nature. When the 
spectrum of the standard star was broadened to the same 
rotational velocity of FR Cnc, the profiles of the CCF be- 
came smoother and could be fitted with a Gaussian profile, 
(see Fig. [9l bottom). 

The irregular profiles of the CCF (double peaks and 
asymmetries) can produce significant errors in RV measure- 
ments (see Fig.|9]). These irregularities may come from pho- 
tospheric activity features on the stellar surface. They can 
distort the profile of the photospheric lines and induce vari- 
ations in the peak of the CCF. However, this behavior may 
be caused by the difference in rotational velo city (v sin i) 
betw een the target and standard star (see e.g. ICalvez et al.l 
l200'iD when the standard is a much slower rotator than the 
target. The CCF is essentially the broadening function that 
would be applied to the template spectrum. Use of a broad- 
ened template removes higher moments introduced from the 
starspots, enabling a relatively unbiased estimate of RV to 
be determined (Fig. E]). A mean RV of 17.8 ± 1.6 is ob- 
tained for the 2004 data set and is in good agreement with 
the RV derived through optimisation of parameters in the 
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Table 7. Stellar parameters of FR Cnc 



T,p 


Vmoan 


B - Vmoan 


^^phot 
(d) 


V sin i 
(km s^l) 


i 
o 


K7V 


10.44 


1.10 


0.826518±0. 000015 


46.2±0.8 


55 



Table 8. Astrometric and kinematic parameters of FR Cnc 



7r 

(marcsGc) 


COS 3 
(marcscc yr~^) 


(marcscc yr~^) 


(km s~l) 


V ±av 
(km s~l) 


(km s~^) 


^'Total 
(km s"^) 


30.24±2.a3t 


-98.1±1.6* 


-91. Oil. 5* 


-19.02±0.62 


-18.99±1.08 


-8.03±1.58 


28.05 



From Hipparcos and Tycho — 2 catalogues |ESA||19971) . 




Figure 9. (Top) An example of CCF of FR Cnc in FOCES04 
observing run. Irregular profiles can be seen in the peak. These 
irregularities can produce significant errors in RV determination. 
(Bottom) The same CCF obtained when we broadened the stan- 
dard star to FR Cnc rotational velocity. Irregular profiles become 
smoother and could be fitted with a Gaussian. 



Doppler imaging process (Section 7), where we find Vhci = 
18.6 ± 0.6. Therefore all the radial velocities given in this 
paper have been calculated by cross-correlation with this 
rotational broadened spectrum of the standard star. 

In Table [9] we list, for each spectrum, the heliocentric 
radial velocities (Vhci) with their corresponding errors (cry) 
obtained as weighted means of individual values deduced 
for each order in th e spectra. We also list data points from 
lUpgren et al.l (|2002l ) for comparison. 

Those orders which contain chromospheric features and 
prominent telluric lines have been excluded when determin- 
ing the mean radial velocity. 



5.3 Rotational Velocity 

By using the program jstarmod (see Sect. 15.1]) we have ob- 
tained the best fits for each observing run using v sin i values 




16 18 20 22 

RV (km/s) 

Figure 10. Bisector velocity span vs. RV for FOCES04 observing 
run. The clear negative correlation indicates that RV variations 
arc due to stellar activity (see Sect. [5T2t . 



of ss35 km s~^. However, this value depends on the rotation 
of the standard star which has non-zero rotation. Therefore 
the obtained value can only be used as an approximation. 

To determine an accurate rotational velocity of 
this star we made use o f the following method (see 
iMartmez-Arnaiz et al.l [20101 for details ). Rotational veloci - 
ties, V s'mi can be written as follows (see lQueloz et all (|l998l ) 
and references therein): 



_ 2 
— <^obs 



2 
^0 



V sini = A\ a 



(3) 



where A is a coupling constant which depends on the spec- 
trograph and its configuration. The spectrum of each of 
these stars was broadened using the program jstarmod 
from V sini = 1 km s~^ up to 50 km s^^ and the respec- 
tive CCF was calculated. A was found for every spectro- 
graph by fitting the relation (usini)^ vs o"oi,s. It is well 
known that (jq is a function of the broadening mechanisms 
which are present in the atmosphere of the star, except ro- 
tation (|Melo et al.ll200j ) . Since the broadening mechanisms 
are function of the temperature and gravity, we may expect 
a dependence of ao on the temperature. To determine this 
dependence we use synthetic spectra with no rotational ve - 
locity computed using the ATLAS9 code bv (lKuruc3 19931) 
adapt ed to work under a Lin ux platform bv lSbordone et al] 
(|2004 ) and ISbordond (|2005h . Once A is determined and ctq 
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Table 9. Radial and Rotational Velocities 



Run 


HJD 

(2400000 +) 


Vhcl 
(km 




(TV 


V SXTLZ 

(km s^l) 


FOCES04 


53098 


3713 


16 


95 


± 


0.58 


44.1 ± 1.9 


FOCES04 


53099 


3411 


18 


42 


± 


0.57 




FOCES04 


53099 


4230 


17 


34 


± 


0.66 




FOCES04 


53099 


4546 


16 


74 


± 


0.69 




FOCES04 


53100 


3156 


16 


62 


± 


0.64 




FOCES04 


53100 


3917 


15 


36 


± 


0.56 




FOCES04 


53101 


3264 


18 


70 


± 


0.48 




FOCES04 


53101 


4495 


20 


97 


± 


0.55 




FOCES04 


53102 


3306 


19 


30 


± 


0.55 




FOCES04 


53102 


4753 


17 


68 


± 


0.55 




FOCES06 


54086 


6770 


19 


85 


± 


0.63 


37.2 ± - 


FOCES06 


54088 


5918 


19 


62 


± 


0.45 




FOCES06 


54091 


6192 


20 


45 


± 


0.53 




FOCES07a 


54156 


4482 


19 


98 


± 


0.52 


41.3 ± 2.7 


FOCESOTa 


54158 


5625 


20 


41 


± 


0.80 




FOCES07b 


54228 


3461 










43.6 ± 3.2 


FOCES07b 


54229 


3496 












FIES08 


54547 


4408 


17 


82 


± 


0.74 




Upgren02 


51626 


695 


27 ± 2.3 




Upgrcn02 


51626 


730 


24 ± 4.1 





calibrated with the color index [B — V), acts (width of the 
CCF of the star when is correlated with itself) is measured 
for each star, usini can be directly calculated using the above 
formula ((3|. 

In Table [2] we list, for each observing run, the averaged 
II sin i value obtained. From Table [9] we estimate uncertain- 
ties of 1.9-3.2 km s~^ based on the standard deviations as 
each measurement epoch. It is likely that the wsini values 
vary by more than these uncertainties since PR Cnc is very 
active and exhibits starspots (see Section 7) that signifi- 
cantly distort the rotationally broadened absorption lines. 
Moreover, while the observations taken in the year 2004 
cover a complete rotation cycle, those usini measurements 
at other epochs, only include one to two observations and 
are likely to yield more biased results (i.e. depending on the 
location of starspots at the observation phases). In Section 
7, we model the line profile using our Doppler imaging code, 
allowing for the presence of spots. The resulting fits are thus 
likely to give a more accurate rotation velocity, wsini, for 
FR Cnc. 



but the Eggen kinematic criteria (see iMontes et al.l (|200ll ) 
for details) are negative, showing that FR Cnc could be not 
a member of any moving group (MG hereafter) in the young 
disc area. 

The classical view of MGs (e.g. Eggen 1984a), i.e. they 
come from the remnant of a star-forming cloud has been 
discussed in recent years. Several studies (e.g. Famaey et al. 
2007, 2008; Antoja et al. 2008; Zhao et al. 2009) seem to 
support a dynamic or resonant mechanism origin. While 
both theories are feasible, we will just take into account 
that the percentage of contamination of the young disc 
space velocity area by ol d field population is high (see 
iLopez-Santiago et al.|[2009l and reference therein) and so age 
constraints are needed to assess if FR Cnc belongs to any of 
young disc moving groups. 

5.5 The Li I A6707.8 line 

Li I A6707.8 spectroscopic feature is an important diagnostic 
tool for assessing the age in late-type stars, since Lithium 
is destroyed easily by thermonuclear reactions in the stellar 
interior. 

The spectral region of the resonance doublet of Li I at 
A6708 A is covered in all the high resolution observations. 
We measured the equivalent width {EW hereafter) in seven 
spectra. 

Due to the small value, we were not able to measure 
EWs directly, or use spectral subtraction technique. There- 
fore EWs have been obtained using the IRAF task SBANDS, 
performing an integration within a band of 1.6 A centred 
in the lithium line (Maldonado et al. 2010). We have ob- 
tained an averaged value of 54 mA. In our high resolution 
spectra the Li I line is blended with Fe I 6707.4 A. To sub- 
tract the_^]e I contribution we used the color-index relation 
from iFavatal (|l993l ). Therefore, we obtain the final value of 
EW(Iji l) to be 34 mA. In Fig. [11] we plot as an example, a 
spectra formed by co-adding the FOCES04 run spectra and 
we indicated the position of the Li I line. 

By comparing this value with stars of similar spectral 
type (K5-K7 with V-I = 1.3) in other MGs members or 
cl usters of a well-known age ( in the same way as in Figure 3 
of lLopez-Santiago et al.|[200£ ). it is in agreement with being 
a young object between 10-120 and compatible with being 
IC 2391 MG member. 



5.4 Kinematics 

We computed the galactic space-velocity components {U, 
V, W) and their ass ociated errors of FR Cnc usin g the pro- 
cedure described by Ijohnson fc SoderblornI (| 19871 ') modified 
by Lopez-Santiago fsee lMontes et al.ll200ll ). This procedure 
uses J2000 coordinates and takes into account correlation in 
the measures of Hipparcos. We use averaged RV calculated 
here (18.63 ± 0.14 km s~^) and the proper motions and 
parallax from Hipparcos. 

The obtained values of the components with its module 
Vrotai and associated errors are given in Table [H] The ve- 
locity components in the {U, V) diagram are clearly within 
the young disc population boundaries (Eggen 1984a, b, 1989; 
Monies et al. 2001a,b) indicating that the star belongs to 
the young disc and that it might also belong to the IC 2391 
moving group, mentioned previously by Pandey et al. (2005), 



5.6 Other age indicators 

Ideally, further constraints on FR Cue's age should be 
calculated to confirm its youth. Popular age approx- 
imations are for example th e X-ray fiux-age relation 
(|Mamaiek fc Hillenbrandl l2008l : equation A3) or the com- 
monly used relation between age and Rhk index, that mea- 
sures chromospheric emission in the cores of the broad chro- 
mospheric Ca ll H fc K lines (see e.g. Noves et al.l Il984l : 
iBaliunas et al. 1 1 19961 : iMamaiek fc Hillenbrandl |2008| ). The 
former relation is valid for stars with spectral types earlier 
than that of FR Cnc and so can not be used. Also, using the 
latter relation, from our Ca ll H fc K fiuxes (see Sect. 16.31 
and Table we obtain a logR^f^ = -3.48 which is outside 
the validity range of the activity-age relation but in this case 
it is compatible with the young age. 
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Figure 11. Spectrum in the Li I line region, resulting of co-adding 
the 10 observed spectra of FOCES04 run in order to increase the 
S/N. The position of the Li I + Fe I and Ca I are marked. 
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6 CHROMOSPHERIC ACTIVITY 
INDICATORS 

Both echelle and long slit spectra analysed in this work 
allowed us to study the behavior of the different indica- 
tors from the Ca ii H & K to the Ca ii IRT lines, which 
are formed at different atmospheric altitudes. The chromo- 
spheric contribution to these features was determined by 
usi ng the spectral subt racti on technique describ ed in detail 
bv lMontes et all (|2000l ') and lGalvez et al.l (|2002l '). 

The excess emission EW of different spectral features 
were measured in the subtracted spectra. In Table [10] we 
give the EW for the Ca ii H & K, He, US, H7, H/3, Ha, 
and Ca 11 IRT (AA8498, 8542, 8662 A) lines for the echelle 
spectra. These EWs were converted to an absolute surface 
fluxes by u sing the empirical stellar flux scales calibrated by 
iHalll I 19961 ) as a function of the star color index. In our case, 
we used the B — V index and the corresponding coefficients 
for Ca II H & K, Ha and Ca 11 IRT, using for He the same 
coefficients as for Ca 11 H & K, and derived the H5, H7 
and H/3 coefficients of flux by carrying out an interpolation 
between the values of Ca 11 H & K and Ha. The logarithm of 
the obtained absolute flux at the stellar surface (logics) in 
ergs cm~^ s~^ for the different chromospheric activity 
indicators is given in Table [TT] 

Fig. [T5] shows representative observations in the Ha, 
and Ca 11 IRT AA8498, 8542 line regions for high resolution 
spectra. Fig. [13] shows representative observations in the Ha 
for low resolution spectra. Fig. [TJ] shows a closer view of one 
spectrum from Figs. [12] and [13] where the emission can be 
better seen. 



Figure 12. Spectra in the Ha (left side) and Ca 11 IRT AA8498, 
8542 (right side) line regions for FOCES04 observing run; clear 
wide and prominent emission arises over the continuum from Ha 
and a clear emission in the core of the absorption line is seen in 
the Ca II IRT lines. 



Ha 



Observed 



HFOSC 2004 




6400 6500 6600 6700 

Wavelength (A) 

Figure 13. Sample of spectra in the Ha line region for HFOSC04 
observing run; clear emission arises over the continuum from Ha 
in this low resolution spectra. 
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Figure 14. A representative spectra of Figs. [121 and 1131 shovfing 
a zoom of tiie Ha (left side) and Ca II IRT (center) line regions 
for FOCES04 and Ha (right) line region for HFOSC04. 
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Figure 15. Variation of Ha EW vs phase in HFOSC04 run. 
Different symbols represent different nights. The second night is 
in triangles, showing higher values of the Ha EWs. 
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Figure 16. Top: Variation of Ha £'14^ vs phase in FOCES04 run. 
Bottom: Variation of Ca II IRT EW vs phase. 

observations. This can be interpreted as an activity cycle 
of at least 4-5 years, similar to Sun or other stars activity 
cycles, but further follow up is needed to confirm this. Ac- 
curate stellar activity cycle can prove useful for study the 
dynamo interface and activity cycle-rotation-spectral type 
mechanisms in the stars (see e.g. Lorente & Montesinos 
2005). 

The persistence of Ha emission indicates that it is a very 
active BY Dra system, but the vast range of variability levels 
make this star unusual and interesting for further study. 



6.1 The Ha line 

We analysed the Ha line region for all the spectra. This line 
in the obtained spectra is always observed in emission above 
the continuum (see Figs. [121 and 113p . 

Measuring the EW of this line, we found that the EW 
average of the Ha emission is quite different in every sea- 
son, showing significant variability in time-scales of a year. 
EW{Ra) = 3.23 A for FOCES04 run while, EW{Ra) = 
1.67, 1.72 and 1.87 A for FOCES06, FOCES07a and FO- 
CES07b respectively, £W(Ha) = 1.72 for the only value of 
FIES08 run. For the low resolution spectra we have an av- 
erage value (in 37 spectra taken during three consecutive 
nights) of EWiRa) = 2.14 A in HFOSC04 run. 

Fig. [m shows the variation of EW vs phase (calcu- 
lated with the photometric period) in HFOSC04 run. Dif- 
ferent symbols represent different nights. The second night 
(triangles) shows higher values of the Ha EWs. Top of Fig. 
1161 represent the variation of EW vs phase in the FOCES04 
run. 

Comparing the variations of Ha EW between the runs 
and the variation in each run when possible, we note ac- 
tivity level variations on a month-long timescale, from one 
year to the next (see Table UOp . Making the comparison 
with the photometry, we notice a correspondence between 
marked variations in the light curve and i?VK(Ha). In the 
year 2004 both photometry and spectroscopy show a high 
level of FR Cnc activity, while it is decreasing rapidly in the 
year 2005 and then remains on that level during our further 



6.2 The H/J, H7 and US lines 

We can see the absorption of H/3, H7 and US Balmer lines 
filled in with emission in the observed spectra. Fig. [17] plots 
a representative subtracted spectra of these three lines in 
dififerent nights of FOCES04 run. 

The variation of these lines with rotational phase and 
from season to season follows the same trend as Ha varia- 
tion. 

We also measured the ratio of excess emission in the 
Ha and H/3 lines ( g^j^^j ) and the ratio of excess emission 
with the correction: 

given bv lHall fc RamsevI (| 19921 ). This corrects the absolute 
flux density in these lines for the color difference in the com- 
ponents. We have obtained a mean value of ^ 2.5. This 
value is in the limit between the presence of prominence-like 
material on or above the stellar surface (jBuzasi, ,19891 and 
iHall fc Ramsevlll992l ). 



6.3 Ca II H & K and He 

The Ca 11 H & K line region is included in most of the spectra 
but the efficiency of the spectrograph and the CCD decreases 
very rapidly due to the position of these lines at the end of 



(4) 
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7 DOPPLER IMAGING 
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Figure 17. Representative subtracted spectra of H/3, H7 and H<5 
lines in different nights of FOCES04 run. 



the echellogram. Therefore, the obtained S/N ratio is very 
low, and the normahzation of the spectra is very difficult. 
In many cases we could not measure the EW Hues and in 
other we measured them in the observed spectra as it was 
not possible to apply the spectral subtraction for this region 
(see Tables [TOl and [TT|). 

Strong emission in the Ca ll H & K is seen despite the 
low S/N. 



6.4 Ca II IRT lines (AA8498, 8542, and 8662) 

The three lines of the Ca ll (IRT) are included in all our 
echelle spectra. In all of them a clear emission arising in the 
core of the absorption lines is seen (see Fig. I12|) . 

Averaged values of emission EW are EW{C& ll) = 
0.74, 0.93 and 0.81 A for AA8498, 8542, and 8662 in FO- 
CES04 run. EW{C& ll) = 0.49, 0.72 and 0.61 A in FO- 
CES06 run. Ewlca. ii) = 0.66, 0.86 and 0.70 A in FO- 
CES07a. EW{C& ii) = 0.79, 1.03 and 0.77 A in FOCES07b 
and EW{Ca. ll) = 0.71, 0.89 and 0.79 A in FIES08. 

The variation of the emission in these lines is significant 
although not as strong as in Ha (see Table (TOjl. 

Only for the FOCES04-run we have enough data points 
to compare variations between Ha EWs. and Ca ll EWs. 
In other active stars a clear anticorrelation is usually seen 
(as their emissions come from different features in th e 
stellar surface ) , see e. g. stars in lArevalo fc Lazarol (Il999h. 



iMontcs et al] (|2000l '). ICalvez et all (|2002h , ICalvez et al] 
(2009) etc., but it is very weak in the case of FR Cnc (see 
Fig. [ISI). 

In addition, we have calculated the ratio of excess emis- 
sion EW , -f^J^, which is also an indicator of the type of 
chromospheric structure, which produces the observed emis- 

1.5-3 are measured, 
9j_the limit of 



£^8498 



sion. In solar plage values of 
while in solar prominence the val ues are 
an optically thin emitting plasma ()Chestej|l99ll l. We have 
found a lafiiiRJl.S, indicating that Ca ll IRT emission comes 

-£'8498 ° 

from plage-like regions. 



Since FR Cnc is a rapid rotator with considerable broad- 
ening of spectral lines, we generated an indirect starspot 
map usin g the Doppler T omography of Stars (DoTS) imag- 
ing code (|Cameronl 1200 ll ) . In order to detect the line dis- 
tortion due to starspots in the high resolu tion spectra, 
we have applied least-sq uares deconvolution (|Donati et al.l 
ll997l : lBarnes et al.|[l998l ') to the 4362-6845 A wavelength re- 
gion of FOCES04 spectra (obtained at 2004 March - April). 
A single line, free of the effects of rotational line blending and 
with high S/N is thus derived. Deconvolution is carried out 
using a T = 2450 K model line list (VALD; see lKupka et al.1 
1999ifc [2OO0l ) which indicates that there are 8345 lines in the 
selected wavelength region with normalized depths of 0. OS- 
LO. Regions around hydrogen Balmer lines, the Mg triplet 
and Na doublet are excluded from the deconvolution. A sin- 
gle line profile with a mean of A = 5460.4 A is derived for 
each observed spectrum. The mean S/N of the input spec- 
tra over the entire 4362A - 6845A wavelength region was 
21.8, while the mean deconvolved line profiles possess S/N 
= 1008 (indicating a gain of 46.2). It should be noted that 
while the profile shown in Fig. [5] (upper panel) represents 
the broadening function of FR Cnc (i.e. rotational veloc- 
ity plus starspot distortions), the LSD profiles in Section 
7 are deconvolved using a linelist rather than a template. 
As such, the LSD profiles still contain the intrinsic stellar 
profile, and any distortions due to starspots will appear to 
possess a lower amplitude when compared with a template- 
derived CCF. 

For the imaging procedure, we used the standard star, 
HD 151877, to represent the local intensity profile of a slowly 
rotating star. A two temperature model with Tphot = 4250 
K and Tspot = 3000 K was used. The starspot image there- 
fore represents the spot filling fact or. Details of the Doppler 
imaging technique can be found in lCameronI (|200lf ). We op- 
timized the goodness of fit to the 10 deconvolved profiles 
for heliocentric RV, axial inclination, equivalent width and 
Dsini, finding Vhei = 18.6 ± 0.6, vsini = 46.2 ± 0.8 km s~^ 
and i = 55 ± 5°. The vsini value is higher than, but still 
consistent with, our mean usini value derived in Section 
5.3. Since we here use all rotation phases to derive our best 
fit Doppler image, it is likely that the result is less biased 
than the previous results in Section 5.3 that do not take 
account of the presence of starspots, which affect the pro- 
file shape. Starspot distortions in the wings of the profiles 
for instance, are likely to lead to underestimations of vsini 
when using a single spectrum, leading to a systematically 
lower mean estimation of vsini. We note however that our 
result for the 2004 observations is consistent with the value 
tabulated in Table 9 (i.e. 44.1 ± 1.9 km s"^). Fig. [TSl shows 
the deconvolved profiles (and phases of observation) while 
Fig. [19] shows mercator projections of the starspot image of 
FR Cnc. 

The surface map in Fig. [19] (upper panel) indicates that 
FR Cnc possesses considerable spot coverage as suggested 
by the time varying distortions in the deconvolved profiles 
in Fig. 1 181 (left panel - spectroscopic data only). For compar- 
ison with the 2003-2004 (JD = 2452934-2453146 in Table 
1) lightcurve presented in Fig. 3 (top right), we have gen- 
erated a synthetic lightcurve based on the spectroscopically 
derived image in Fig. [19] (upper panel). While we are able 
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to recover the observed Ughtcurve morphology, we consid- 
erably underestimate the amplitude (Fig. 1201 upper panel). 
This is likely due, at least in part, to loss of information in 
the Doppler imaging process (i.e. of the spectroscopic data) 
owing to resolution and noise constraints imposed by the 
fitting procedure which must optimise both the fit and the 
image entropy. With higher S/N ratio data and spectral res- 
olution it is probable that more precise fitting would yield 
spot filling factors indicating cooler spots. In addition, there 
is a tendency for spots close to the equator (<20°) to be re- 
constructed at slightly higher latitudes than they are phys- 
ically situated. This is due to the smaller entropy penalty 
of placing a smaller spot away from the equator as com- 
pared with a larger spot nearer the equator. The resulting 
spot distribution will naturally lead to predicted lightcurves 
with underestimated amplitudes since equatorial spots yield 
higher amplitudes. 

We find a spot filling factor of 6 per cent from our spec- 
troscopic Doppler image. This is at the lower end of the typ- 
ical values (up to ~ 30 per cent) found for active stars by 
O'Neal, Saar & Neff (1996), O'Neal, Neff & Saar (1998) and 
O'Neal et al. (2004) using TiO indicators. We have therefore 
reconstructed an image that also makes use of the 2003/2004 
y-band photometric data (e.g. see Barnes et al. 1998). It 
should be noted that while the lightcurve data spans sev- 
eral months, the spectroscopic data were taken over only a 
few days. Since the photometric observations taken within 
± 1 month of the spectroscopic data do not show signif- 
icant deviation from the photometric dataset as a whole, 
we made use of all photometric points. The resulting simul- 
taneous fits to the data are shown in Fig. [18] (right panel) 
and Fig. [20] (bottom panel), while the reconstructed image is 
shown in Fig. [T5] (bottom panel). Simultaneous imaging that 
makes use of spectroscopic and photometric data must strike 
a balance, such that reasonable fits to both data sets are 
achieved. We have already outlined (above) possible modest 
biases in an image derived from a spectroscopic fit to finite 
S/N. Lightcurve data alone does not contain sufficient infor- 
mation to make reliable images of single stars. Any complex 
spot structure or spot groups will be rendered simply as a 
single spot. Photometric imaging favours reconstruction at 
lower latitudes because a smaller spot is more easily able to 
reproduce the amplitude variation (and yields a lower image 
entropy penalty) in an observed lightcurve than a slightly 
larger spot at higher latitude. Allowing the photometric fit 
to more closely match the observed lightcurve results in 
a slightly poorer fit to the spectroscopic data. The fitted 
lightcurve now closely matches the observed lightcurve, since 
more large scale (albeit unresolved) starspot structure has 
been recovered. This has necessitated some loss of image de- 
tails (Fig. 1191 lower panel) as compared to the spectroscopic 
reconstruction alone (Fig. 1191 upper panel). The filling factor 
is now greater at 9 per cent, adding further weight to the 
argument that spectroscopy alone is unable to recover all 
starspot information. Indeed, if the results from TiO are to 
be believed, there may be further unresolved global starspot 
covereage on FR Cnc and all other stars that have been sub- 
jected to such studies. 

We also carried out a bisector analysis on the LSD pro- 
files, calculating the bisectors in the same way as in Sect. 
5.2. As we have mentioned, the LSD profiles should show 
a smoothing of the bisectors, reflected in the lower value of 
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Figure 18. The figure shows the deconvolved profiles, with the 
phases of observation. The variation of the profile shape, due to 
starspots rotating into and out of view, is clearly seen. 



the bisector we obtained comparing with the original spec- 
tra. As expected, the anticorrelation is still seen, but here 
we have a strongest Pearson correlation coefficient of r = 
-0.8992. 

As might be expected for such an active star, FR Cnc 
exhibits a high degree of spot coverage. No spots are visi- 
ble in the phase range 0.1-0.3. This agrees well with the 
photometric l^-band light curve plotted for the same pe- 
riod (Fig. [3] 2003-2004) which shows a maximum at phase 
~ 0.2. This contrast with a significant degree of spotness 
at other longitudes is likely responsible for the high degree 
of photometrical modulation seen in this star. The starspot 
distribution of FR Cnc is also of interest since it is one of the 
latest spectral types to have bee n imaged. Comparis on with 
the K5V rapid rotator LO Peg, [Barnes et al.l (|2005l '). shows 
a significantly different spot distribution. While LO Peg ro- 
tates more rapidly (v sin i = 65.9 km s~^), showing predom- 
inantly mid-high latitude spot structures and a polar spot, 
FR Cnc shows only mid-latitude starspots. Small scale spot 
struct ure is expected to b e variable on a day-long timescale, 
while [Jeffers et al] \2m% have shown that the polar spot 
on the KOV star, AB Dor, is variable in extent over periods 
of years. While G, and K stars such as AB Dor often show 
strong polar caps, no such features were seen on the M1-M 2 
dwarf stars HK Aqr and RE 1816-1-541 ([Barnes et al.ll200ll '). 
This change in starspot location (i.e. with latitude) may be 
due to a change from a convective shell-type dynamo to a 
fully convective dynamo. Predominantly mid-high latitude 
spots are expected to arise from a rapidly rotating s olar-like 
dynamo ([Schiissler et al.lll996l : [Cranzer et al.|[2000l ) while a 
distributed dynamo may be expected to produce spots at all 
latitudes. However, only further spectroscopic time series ob- 
servations of FR Cnc, with a higher cadence of observations 
(enabling better surface resolution to be obtained), would 
reveal the long term stability of the starspot patterns. 



8 SUMMARY AND CONCLUSIONS 

We have carried out a photometric, polarimetric and spec- 
troscopic study of FR Cnc. 

Optical B-band photometry was carried out during 



FR Cnc Revisited: Photometry, Polarimetry and Spectroscopy 15 



Cnc - March/April 2QQ4 (Spect) 



FR Cnc - h 




0,1 0,2 0,3 0,4 0,5 0,6 0,7 0.8 



phase 



FR Cnc - March/April 2004 (Spect) and 2005/2004 (Phot) 



0,1 0,2 



0.4 0.5 0.6 0,7 0,8 0.9 110 




phose 



Effective spot area 



Figure 19. Mcrcator projection of the starspot image of FR Cnc. A high degree of spot coverage is revealed. There are no visible spots 
in the phase range 0.1-0.3. FR Cnc shows only mid-latitude starspots and no polar spot like other typical late K stars. The vertical 
lines mark the corresponding phase for 10 spectra we have used for deconvolution. The right hand plot shows how the mean effective 
spot area changes with latitude. 



Table 10. EW of Chromospheric Indicators 
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^ Values measured with low S/N. 
^ Values not measured due to the low S/N. 
* Values measured in the observed spectra (spectral subtraction is not applied). 



2007 March - 2008 February at Terskol Branch of the As- 
tronomy Institute (Russia), There are no pecuhar features 
in the 2008 hght curve, while in the 2007 photometry two 
brightening episodes were detected. One of them occurred 
at tlie same pliase as the flare of November 23, 2006 (phase 
= 0,88) and probably indicates that both of these events 
(i,e, the flare on November 23, 2006 and the photometric 
brightening episode) originated from the same long-living 
active region on FR Cnc, The non-detection of any other 
flares in our photometry except 2006 November 23 implies 
that FR Cnc has a low frequency of flares. 



We analysed ASAS-3 photometry obtained in 2002- 
2008 in l/-band. No evidence of flares in ASAS-3 data were 
found. The profiles of variability are different from season 
to season. The mean magnitude in V-band remained the 
same (Vmoan = 10,439 mag) during 2002-2008, while the 
amplitude decreased abruptly in 2005, The proposed inter- 
pretation is a redistribution of spots/spot groups over the 
surface of the star, while the total percentage of the spot- 
ted area was assumed to be constant within the error limits, 
A detailed periodogram study of the ASAS-3 photometric 
data enabled us to derive a more accurate value for the pe- 
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Table 11. Emission Fluxes 
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Notes as in previous tabic. 
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Figure 20. Photometric lightcurvcs calculated from the surface 
maps and superimposed on ylSj45-5-lightcurve for corresponding 
season of observations. 



riod of FR Cnc. We find that P = 0.08265 ± O.OOOOf 5 d. In 
addition, we also presented BVRcIc piiotometric calibration 
of 166 stars in FR Cnc vicinity, whose V-magnitude is in the 
range of 9.85-18.06 mag. 



The BVR broad-band polarimetric observations of 
FR Cnc have been obtained at ARIES in Nainital (India) at 
Manora Peak. The observed polarization in B-band is well 
matched with the theoretical values expected for Zeeman 
polarization model. However, the observed polarization in 
V and R bands slightly exceeds the theoretical values and 
Thompson and Rayleigh scattering from inhomogeneous re- 
gions are not enough to explain the observed polarization 
excess. Therefore the excess of linear polarization should 
come from an additional source of polarization. Taking into 
account that we conclude that FR Cnc is not a binary, the 
mechanism which can produce additional linear polarization 
is probably scattering in circumstell ar material d istributed 
in an asymmetric geometry (e.g. see IPandev et al.. . i2009. 'l . 

A total of 58 spectra of FR Cnc, which have been ob- 
tained in 2004-2008, were analysed in this work. Based on 
our spectroscopic observations, FR Cnc was classified as 
K7V star. RV analysis supports the single nature of FR Cnc. 
Anticorrelation between BIS and RV also indicates that the 
RV variations are due to stellar activity variations and not 
due to a secondary companion. 

The kinematics study, based on obtained galactic space- 
velocity components {U, V, W) of FR Cnc, shows that this 
star clearly lies in the young disc population velocity space 
and might also belong to IC 2391 moving group, although 
the Eggen kinematic criteria shows that FR Cnc may not be 
a member of any MG in the young disc area. The Li I A6707.8 
averaged EW measured is 34 mA, giving the spectral type 
of FR Cnc, it is in agreement with being a young object 
between 10-120 Myr. 

The Ha line was always observed above the continuum 
in all the obtained spectra. Measuring the EW of this line, 
we found that the Ha emission EW average in every season 
is quite different. In 2004, as with the photometry, spec- 
troscopic indicators of chromospheric activity show a high 
level of activity which decreased in 2005. The Ca 11 (IRT) 
is included in our echelle spectra. From the ratio of excess 
emission EW we found that in FR Cnc, Ca 11 emission comes 
from plage-like regions. We noticed that FR Cnc can show 
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an activity cycle of 4-5 years, although further follow up will 
confirm this periodicity. 

Since FR Cnc is a rapid rotator, we generated an indi- 
rect starspot map using the Doppler Tomography of Stars 
imaging code. From it we derive v sin i = 46.2 ± 0.8 km s^^ 
and i = 55 ± 5° . FR Cnc belongs to one of the latest spectral 
types to have been imaged with the Doppler Tomography. 

We independently estimated a rotational velocity of 
FR Cnc during our observations using Queloz method (Sect. 
5.3) and by the D.I. fits (Sect. 7). Although they are differ- 
ences in the results of the two methods, they are consistent. 
In TableQwe have put as the rotational velocity of FR Cnc, 
the value obtained in the D.I. as it is likely more accurate. 

Despite the short rotation period and its late spectral 
type, FR Cnc shows very few flare events. It shows high 
level of activity as it is a young star, but an unusually short 
variability due to the redistribution of activity features on 
the stellar surface. While this variability is reflected in the 
changes of the amplitude of brightness, the mean brightness 
permanently is nearly constant, indicating that the percent- 
age of stellar surface covered by spots is also constant. The 
spots location is also unusual, not showing a polar spot like 
other F-K stars do but a distribution more resembling those 
seen in M1-M2 dwarfs. Although this may be indicative of a 
distributed dynamo, the mid-high latitude spot locations are 
more suggestive of an interface dynamo under the action of 
rapid rotation. We can only speculate as to whether FR Cnc 
is representative of a regime in which a convective-shell-type 
dynamo gives way to a fully convective dynamo. Polarimet - 
ric observations of the m agnetic field bv lDonatT et all (|2008l ') 
and lMorin et al.l (|2008l ) for example suggest that this occurs 
at a later spectral type of M4, whereas other chromospheric 
indicators s how no obvious change s until later M spectral 
types (e.g. iMohantv fc Basril |2003| ). Further spectroscopy 
with a higher cadence would enable more detailed maps to 
be derived, with multiple epochs enabling the evolution of 
starspots to be investigated. 



ACKNOWLEDGMENTS 

A. Golovin is thankful to Dr. Anju Mukadam and Dr. 
Paula Szkody for useful discussions concerning Fisher Ran- 
domization Test in periodogram analysis, to Dr. Andras 
HoU (Konkoly Observatory, Budapest, Hungary) for valu- 
able comments on VOTable format and help with converting 
photometric sequence to it in order to align on DSS image 
in ALADIN, to Nick Malygin and Dr. Ludmila Pakuliak for 
useful discussions and proof-reading of the manuscript. 

. M.C.Galvez-Ortiz acknowledges the financial sup- 
port from the European Commission in the form 
of Marie Curie Intra European Fellowship (PIEF-GA- 
2008-220679) and the partial support by the Span- 
ish MICINN under the Consolider-Ingenio 2010 Pro- 
gram grant CSD2006-00070: First Scienc e with the GTC 
I http://www.iac.es/consolider-ingenio-gtc ). M.C. Galvez- 
Ortiz and J. Barnes also has received support from RoPACS 
during this research, a Marie Curie Initial Training Network 
funded by the European Commission's Seventh Framework 
Programme. This work was partly supported by the Span- 
ish Ministerio de Ciencia e Innovacion (MICINN), Programa 
National de Astronomia y Astrofisica under grant AYA2008- 



00695, and under grant AYA2008-06423-C03-03, and the Co- 
munidad de Madrid under PRICIT project S2009/ESP-1496 
(AstroMadrid). A. Golovin is thankful to Universidad Com- 
plutense de Madrid for hospitality and for all the efforts 
and the help during his visit to Spain in 2008 July and 2009 
February. 

This publication made use of the ALADI N interactive sky 
atlas , operated at CDS, Strasbourg, France jBonnarel et al.l 
l2000t ) and of NASA's Astrophysics Data System. 



REFERENCES 

Alekseev I. Y., 2003, ARep, 47, 430 

Antoja T., Figueras F., Fernandez D., Torra J., 2008, A&A, 
490, 135 

Arevalo M. J. & Lazaro C, 1999, AJ, 118, 1015 
Bahunas S. L., Nesme-Ribes E., Sokoloff D., Soon W. H., 

1996, ApJ, 460, 848 
Barden S. C, 1985, ApJ, 295, 162 

Barnes J. R., Collier Cameron A., Unruh Y. C, Do- 

nati J. F., Hussain G. A. J., 1998, MNRAS, 299, 904 
Barnes J. R., Collier Cameron A., 2001, MNRAS, 326, 950 
Barnes J. R., Cameron A. C, Lister T. A., Pointer G. R., 

Still M. D., 2005, MNRAS, 356, 1501 
Beavers W. I., Fitter J. J., Ketelsen D. A., Oesper D. A., 

1979, PASP, 91, 698 
Bessell M. S., 1979, PASP, 91, 589 
Bonfils X. et al., 2007, A&A, 474, 293 
Bonnarel F. et al., 2000, A&AS, 143, 33 
Breger M. et al., 1993, A&A, 271, 482 
Breger M. et al., 1999, A&A, 349, 225 
Buzasi D. L. 1989, PhD Thesis, Pennsylvania State Univ. 
Collier Cameron A., 2001, 'Astrotomography - Indirect 

Imaging Methods in Observational Astronomy', Springer 

(Lecture Notes in Physics), 183 
Chester M. M. 1991, PhD Thesis, Pennsylvania State Univ. 
Donati J.-F., Semel M., Carter B., Rees D. E., Collier 

Cameron A., 1997, MNRAS, 291, 658 
Donati J.-F. et al., 2008, MNRAS, 390, 545 
Dorren J. D., Guinan E. F., Dewarf L. E., 1994, ASPC, 64, 

399 

Eggen O. J., 1984a, AJ, 89, 1358 
Eggen O. J., 1984b, ApJS, 55, 597 
Eggen O. J., 1989, PASP, 101, 366 
Elias N. M., II, Dorren J. D., 1990, AJ, 100, 818 
Elvis M., Plummer D., Schachter J., Fabbiano G., 1992, 
ApJS, 80, 257 

ESA, 1997, The Hipparcos and Tycho Catalogues, ESA 
SP-1200 

Famaey B., Pont F., Luri X., Udry S., Mayor M., Joris- 

sen A., 2007, A&A, 461, 957 
Famaey B., Siebert A., Jorissen A., 2008, A&A, 483, 453 
Favata F., Barbera M., Micela G., Sciortino S., 1993, A&A, 

277, 428 

Fisher R.A., 1935, Design of Experiments, Oliver and Boyd, 
Edinburgh. 

Fisher R.A., 1936, J.R.Anthropol.Inst., 66, 57. 
Galvez M. C. et al., 2002, A&A, 389, 524 
Galvez M. C, 2005, PhD Thesis, Universidad Complutense 
de Madrid 

Galvez M. C. et al., 2007, A&A, 472, 587 



18 A. Golovin et al. 



Galvez M. C. et al., 2009, AJ, 137, 3965 

Golovin A., Pavlenko E., Kuznyetsova Y., Krushevska V., 

2007, IBVS, 5748, 1 
Granzer T., Schiissler M., Caligari P., Strassmeier K. G, 

2000, A&A, 355, 1087 
Hall J. C, 1996, PASP, 108, 313 
Hall J. C, Ramsey L. W., 1992, AJ, 104, 1942 
Huovelin J., Saar S. H., 1991, ApJ, 374, 319 
Jeffers S. V., Donati J.-F., Collier Cameron A., MNRAS, 

375, 567 

Johnson H. L., 1966, ARA&A, 4, 193 
Johnson D. R. H., Soderblom D. R., 1987, AJ, 93, 864 
Kallinger T., Reegen P., Weiss W. W., 2008, A&A, 481, 
571 

Kazarovets E. V., Samus N. N., Durlevich O. V., 
Frolov M. S., Antipin S. V., Kireeva N. N., Pas- 
tukhova E. N., 1999, IBVS, 4659, 1 

Kozhevnikova A. V., Alekseev I. Y., Heckert P. A., 
Kozhevnikov V. P., 2006, IBVS, 5723, 1 

KupkaF., PiskunovN., Ryabchikova T. A., Stempels H. C, 
Weiss W. W, 1999, A&AS, 138, 119 

Kupka F. G., Ryabchikova T. A., Piskunov N. E., Stem- 
pels H. C, Weiss W. W, 2000, BaltA, 9, 590 

Kurucz R. L. 1993, in Astronomical Society of the Pacific 
Conference Series, Vol. 44, lAU Colloq. 138: Peculiar ver- 
sus Normal Phenomena in A-type and Related Stars, ed. 
Dworetsky M. M., Castelli P., Faraggiana R., 87. 

Kuschnig R., Weiss W. W., Gruber R., Bely P. Y., 
Jenkner H., 1997, A&A, 328, 544 

Lenz P., Breger M., 2005, Communications in Asteroseis- 
mology, 146, 53 

Lopez-Santiago J., Monies D., Fernandez-Figueroa M. J., 
Ramsey L. W., 2003, A&A, 411, 489 

Lopez-Santiago J., Micela G., Monies D., 2009, A&A, 499, 
129 

Lopez-Santiago J., Monies D., Galvez-Ortiz M. C, 
Crespo-Chacon I., Martmez-Arnaiz R. M., Fernandez- 
Figueroa M. J., de Castro E., Cornide M., 2010, A&A, 
514, A97 

Lorente R., Montesinos B., 2005, AJ, 632, 1104 

Maldonado J., Martmez-Arnaiz R. M., Eiroa C, Monies 
D., Montesinos B., 2010, A&A, 521A, 12 

Mamajek E. E., Hillenbrand L. A., 2008, ApJ, 687, 1264 

Martmez-Arnaiz R., Maldonado J., Monies D., Eiroa C, 
Montesinos B., 2010, A&A, 520, A79 

Martinez Fiorenzano A. F., Gratton R. G., Desidera S., 
Cosentino R., Endl M., 2005, A&A, 442, 775 

Medhi B. J., Maheswar G., Brijesh K., Pandey J. C, Ku- 
mar T. S., Sagar R., 2007, MNRAS, 378, 881 

Melo C, Pasquini L., de Medeiros J. R. 2004, in lAU Sym- 
posium, Vol. 215, Stellar Rotation, ed. Maeder A., Ee- 
nens P., 455 

Mohanty S., Basri G., 2003, ApJ, 583, 451 

Monies D., Fernandez-Figueroa M. J., De Castro E. et al., 
2000, A&AS, 146, 103 

Monies D., Lopez-Santiago J., Galvez M. C. et al., 2001a, 
MNRAS, 328, 45 

Monies D., Lopez-Santiago J., Fernandez-Figueroa M. J., 
Galvez M. C, 2001b, A&A, 379, 976 

Montgomery M. H., Odonoghue D., 1999, DSSN, 13, 28 

Morin J. et al., 2008, MNRAS, 390, 567 

Noyes R. W., Hartmann L. W., Baliunas S. L., Dun- 



can D. K., Vaughan A. H., 1984, ApJ, 279, 763 
O'Neal D., Saar S. H., Neff J. E., 1996, ApJ, 463, 766 
O'Neal, D., Neff, J. E. & Saar, S. H. 1998, ApJ, 507, 919 
O'Neal, D., Neff, J. E., Saar, S. H., & Cuntz, M. 2004, AJ, 

128, 1802 

Pandey J. C, Singh K. P., Sagar R., Drake S. A., 2002, 

IBVS, 5351, 1 
Pandey J. C, 2003, BASI, 31, 329 

Pandey J. C, Singh K. P., Drake S. A., Sagar R., 2005, 
AJ, 130, 1231 

Pandey J. C, Medhi B. J., Sagar R., Pandey A. K., 2009, 

MNRAS, 396, 1004 
Panov K., Goranova Y., Genkov V., 2000, IBVS, 4917, 1 
Pepper J., Stanek K. Z., Pogge R. W., Latham D. W., 

DePoy D. L., Siverd R., Poindexter S., Sivakoff G. R., 

2008, AJ, 135, 907 
Ferryman M. A. C. et al., 1997, A&A, 323, L49 
Pojmanski G., 2002, AcA, 52, 397 

Queloz D., AUain S., Mermilliod J.-C, Bouvier J., 

Mayor M. 1998, A&A, 335, 183 
Queloz D. et al., 2001, A&A, 379, 279 
Saar S. H., Huovelin J., 1993, ApJ, 404, 739 
Sbordone L., Bonifacio P., Castelli F., Kurucz R. L., 2004, 

MSAIS, 5, 93 

Sbordone L., 2005, Memorie della Societa Astronomica 

Italiana Supplement, 8, 61 
Schachter J. F., Remillard R., Saar S. H., Favata F., 

Sciortino S., Barbera M., 1996, ApJ, 463, 747 
Schmidt G. D., Elston R., Lupie O. L., 1992, AJ, 104, 1563 
Schiissler M., Caligari P., Ferriz-Mas A., Solanki S. K., 

Stix M., 1996, A&A, 314, 503 
Tonry J., Davis M., 1979, AJ, 84, 1511 
Upgren A. R., Sperauskas J., Boyle R. P., 2002, BaltA, 11, 

91 

Voges W. et al., 1999, A&A, 349, 389 

Zhao J., Zhao G., Chen Y., 2009, ApJ, 692, L113 

This paper has been typeset from a T^jX/ ffl^gjX file prepared 
by the author. 



